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EFFECT OF SHALLOW MERIDIONAL CURVATURE ON THE 


VIBRATION OF NEARLY CYLINDRICAL SHELLS* 

By Paul A. Cooper 
Langley Research Center 

SUMMARY 

The effect of the meridional curvature on the minimum natural frequencies of 
nearly cylindrical shells with either positive or negative Gaussian curvature is inves- 
tigated. The positive Gaussian curvature shells have minimum natural frequencies well 
above those of corresponding cylindrical shells. The minimum natural frequencies of 
the negative Gaussian curvature shells generally are below those of the corresponding 
cylinders and evidence wide variations in value, large reductions in magnitude occurring 
at certain critical values of curvature. The large reductions in magnitude result from a 
condition of essentially inextensional vibration. It is demonstrated that membrane theory 
affords a simple method of determining the modal wavelength ratio and curvature combi- 
nations at which these reductions occur, whereas pure bending theory gives a good esti- 
mate of the magnitudes of the associated frequencies. 

INTRODUCTION 

Shells of double curvature are common elements in aerospace vehicle structures. 
Because of the complexity of the equations of motion governing their dynamic behavior, 
numerical methods are commonly employed for solution. (For example, see refs. 1 
to 3.) Numerical methods are adequate and yield accurate results for many applications; 
however, they are not practical for extensive parameter studies of the behavior of shell 
systems even with present-day high-speed computers. On the other hand, if certain 
approximations are made to the equations of motion, closed-form solutions become pos- 
sible in some instances and facilitate rapid calculations. The purpose of the present 
paper is to develop closed-form solutions of approximate equations for the natural vibra- 
tions of a specific class of shells of double curvature in order to determine the effect of 
the meridional curvature on the natural frequencies of vibration of these shells. 

*The information presented herein is based on part of a thesis entitled "Vibration 
of Stressed Shells of Double Curvature," offered in partial fulfillment of the requirements 
for the degree of Doctor of Philosophy in Engineering Mechanics, Virginia Polytechnic 
Institute, Blacksburg, Virginia, June 1968. 



A set of nonlinear homogeneous equations governing the dynamic behavior of pre- 
stressed doubly curved shells of revolution with positive or negative shallow meridional 
curvature is developed. The procedure used in developing these equations is similar to 
that used in references 4 and 5 where Donnell-type shallow shell equations are derived. 

In the present paper, however, the shallowness assumption is made only in relation to 
the meridional direction. By imposing the assumption of zero static meridional rotation, 
linear vibration equations with constant coefficients are obtained. These approximate 
equations govern the infinitesimal natural vibrations of prestressed nearly cylindrical 
shells of revolution and have closed-form solutions. They are solved here for the natural 
vibration of shells with freely supported edges. 

Results from the approximate analysis are compared with more accurate results 
obtained by the use of techniques and equations given in reference 3 to indicate the range 
of application of the present analysis. The approximate analysis is then used to inves- 
tigate the effects of meridional curvature on the natural vibrations of freely supported 
shells as the length and thickness parameters are varied. 

The equations of the membrane and inextensional shell theories are also solved for 
the natural vibration of freely supported shells. The solutions to the membrane theory 
give a simple relationship for determining vibration modes in which shells of negative 
meridional curvature undergo predominantly inextensional deformations, whereas the 
solutions to the inextensional theory yield good approximations to the frequencies of these 
modes. 
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SYMBOLS 


extensional stiffness (eqs. (A14)) 
central rise of shell meridian 
bending stiffness (eqs. (A14)) 

Young's modulus of elasticity 
shell thickness 

nondimensional meridional curvature, 
number of axial half-waves 
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moment resultants associated with vibration state 
number of circumferential waves 
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stress resultants associated with vibration state 
static stress resultants 
surface loading 

circumferential radius of doubly curved shell 

radius of cylinder; circumferential radius at midlength of doubly curved shell 
meridional radii of curvature 
total meridional arc length 
time 

displacement variables defining the vibration state in meridional (|), cir- 
cumferential ( 6 ), and normal directions, respectively 

coordinate measured normal to middle surface 

ratio of circumferential to axial wavelength 

middle-surface strains associated with vibration state 

circumferential coordinate 

middle-surface bending strains associated with the vibration state 
thickness parameter, h/R 
coefficients of equation (B2) (j = 0, 1, 2) 

Poisson's ratio (taken as / u - 0.3 for all calculations) 
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mass density 
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Notations used to identify load and deformation variables: 

(~) indicates variables associated with deformation measured with respect to 

an unstressed shell of double curvature 

(") indicates variables associated with total deformation measured with respect 

to an unstressed undeformed cylinder 

(-) indicates variables associated with the prestress state only 

( ) Q indicates variables associated with initial displacements 

(_) indicates terms not appearing in Donnell-type theory 

Unmarked variables indicate variables associated with the vibration state only. 

DESCRIPTION OF ANALYSIS 


meridional coordinate 


percent ratio of meridional rise to length, §(100) ~ -f^— (100) 

S 8R.£ 

membrane stress function defined in equations (B5) 


natural circular frequency 

R 2 o>2y(i - u 2) 


frequency parameter, 


E 


biharmonic differential operator, r4 + 2R 2 — | + — — 
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To determine the effect of double curvature on the natural vibration behavior of 
shell structures, the class of shells of revolution shown in figure 1 is investigated. These 
shells, chosen for their simplicity, have a constant meridional curvature R| throughout 
their length. The equation of the circumferential radius r which defines the meridional 
shape is 


r = R + R^ 




( 1 ) 
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where R is the circumferential radius at midlength of the shell, R| is the constant 
meridional radius of curvature, s is the total length of meridian, and £ is the merid- 
ional length coordinate measured with the origin at one edge of the shell. Accurate 
results can be found for the natural frequencies of specific shells described by equa- 
tion (1) by the use of numerical procedures such as that given in reference 3. However, 
even with high-speed computers, it is impractical to attempt an extensive parameter 
study of the shell behavior. To achieve a more rapid and economical analysis of the 
effects of various parameters governing the dynamic behavior of the shell configurations 
described by equation (1), an approximate set of equations of motion is developed in 
appendix A. These approximate equations are based on the assumption that the shells 
are shallow in the meridional direction. In this context, equation (1) is replaced by a 
simpler expression in which only the first two terms of a power series expansion for r 
is retained, that is, 


r = R - 



( 2 ) 


Solutions to the approximate equations are obtained in appendix B. The accuracy of these 
solutions is evaluated by comparison of results with more accurate results found by using 
the numerical procedure of reference 3. When the numerical procedure is employed, the 
meridional shape is described by equation (1) and 200 equally spaced intervals are used 
in a finite-difference scheme along the meridian. 

To develop the approximate governing equations, the nonlinear strain-displacement 
relations given by reference 6 are specialized for a cylinder and modified slightly by 
using a criterion of consistency with first approximation shell theory as developed in 
reference 7. The cylinder is then given a small axisymmetric deformation and the 
strains caused by this initial geometric perturbation are removed to establish an 
unstrained state in the resultant doubly curved configuration. The nonlinear equilibrium 
equations corresponding to this unstrained state are developed and linearized. 


The resulting linearized differential equations are 


(I- 


ixV',s \ 2w ,eee (2 - n)x 
+ R^nr ' TI 2 12 w ,me - -§■ v >tt 


W ttfl- ~V„ = 0 




} (3a) 
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and the resulting boundary conditions at £ = Constant are 


u ,0 . L . x 2 |„ X' 


R + { 1 + T v ,rT w ,er° 


ii[-<2 - 4 )v (e + R 2 " 4{{ + (2 - M>w e8| ] - ^w >{ . 0 

w ,«? + ^,o- v ), e = 0 


or u = 0 


or v = 0 


or w = 0 


or w = 0 

j 


The exact solution of equations (3) for freely supported boundary conditions 
= v = w = = Oj, is given by 


u = U mn cos sin n<9 e 


”\ 
iwt 


v = V mn sin cos n0 e ic °t > 

b ■ 

w = W mn sin m < J ^ sin n0 e^ w ^| 


(3b) 


(4) 


where U mn , V mn , and W mn are constant modal amplitudes, m is the number of 
meridional half-waves, n is the number of circumferential waves, and c o is the natural 
circular frequency. Substitution of equations (4) into equations (3a) leads to the following 
characteristic equation for natural frequencies: 


-O 6 + A 2 fi 4 - AjfT + A 0 = 0 


where S2 is a nondimensional frequency given by 


n = R co^L__ii5 


(5) 


(6) 


The coefficients Aj depend upon the shell geometry and mode shape and are given by 
equations (B3). 

The limiting cases of membrane and pure bending (inextensional) behavior are also 
solved in appendix B for freely supported shells. The membrane frequency, with in-plane 
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inertias neglected and with zero initial static stresses, is given by 


^membrane - ~ ~ ^ 


e + » 2 ) 

(i + ef 
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where k| = — is the ratio of the circumferential radius at midlength to the meridional 

radius of curvature and /3 = is the ratio of the circumferential to axial wavelength 

of the vibration mode. It is obvious from equation (7) that the membrane frequency 
vanishes for 


k 



( 8 ) 


The pure bending frequency, with in-plane inertial forces neglected and with zero 

initial static stresses, is given by 

s V u n ( 1 + e 2 ) 2 ” 4 < 9 > 

v. 

where X = — is the ratio of thickness to circumferential radius. The pure bending fre- 
ft 

quencies are independent of meridional curvature and are equal to the natural frequencies 
of a simply supported rectangular plate with aspect ratio 7 tR/s. 

In order to give a more descriptive indication of the degree of shallowness of the 
meridian than that provided by the ratio of radii of curvature an auxiliary param- 
eter r, the percent ratio of the meridional rise to length, is introduced. The central 
rise of the shell meridian (c in fig. 1) is given by 


c = rJi - cos 


2Ri 


( 10 ) 


The central rise can be approximated by the first term of a series expansion of the right- 
hand side of equation (10) if is large compared with s so that r may be repre- 

sented approximately by 


r = §(100) 


8R£ 


( 100 ) 


(ID 


which is consistent with the accuracy of the equations developed in this paper. 


ACCURACY OF SOLUTIONS 


The accuracy of the approximate theory is established by comparing results obtained 
with this theory with corresponding results obtained with the more accurate numerical 
method of reference 3. This comparison is made for several specific doubly curved 
unstressed shells with the same length but different values of k| and results are shown 
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in figure 2. In figure 2, the lowest frequencies are plotted as a function of successive 
circumferential mode numbers for various positive and negative Gaussian curvature 
shells. The solid symbols denote results of the approximate theory and the open symbols 
denote results from the numerical solution. The cylinder results, found by using the 
approximate theory, are shown by a dashed line and are repeated on each figure for pur- 
poses of comparison. For the particular case of an unstressed cylinder, the approximate 
method yields exact solutions since the only approximation is neglect of small nonlinear 
terms related to initial static stress prior to linearization. 

For the slightly curved shells (r = ±1.87%) of figure 2(a), the approximate theory is 
shown to be very accurate, a larger percentage error being evident for the negative 
Gaussian curvature shell than that for the positive Gaussian curvature shell. The accu- 
racy diminishes as the curvature of the shell meridian increases as is shown in fig- 
ures 2(b) (t = ±3.75%) and 2(c) (r = ±5.62%). The positive curvature shell results never 
differ from the more accurate results by more than 8 percent. The negative curvature 
shell results have larger errors but still indicate the general character of the more accu- 
rate solution. 

The lowest natural frequency for the axisymmetric vibration mode n = 0 is asso- 
ciated with a pure torsional mode and is independent of the meridional curvature in the 
approximate theory. The more accurate results, however, indicate that this frequency 
is dependent on curvature to a small degree as can be seen by inspection of figure 2. 

The approximate procedure has a modal solution of sinusoidal form along the 
meridian (see eqs. (4)) whereas the mode shape in the numerical procedure is calculated 
once the frequency is determined and need not necessarily be sinusoidal. Plots of the 
normal deflection shape w for the minimum frequency determined by using the numeri- 
cal procedure are given in figure 3 for particular values of n for the negative curvature 
shells of figure 2. In each case, the number of axial half-waves m determined by the 
approximate theory agrees with the modal shapes given in figure 3. As the curvature 
increases, however, the mode shape begins to deviate from sinusoidal form; thus, the 
sinusoidal modal solution of the approximate theory becomes a less accurate representa- 
tion of the true modal configuration. 

A comparison of the results found by the approximate theory and the numerical 
procedure for a long negative Gaussian curvature shell with a rise- length ratio greater 
than 6 percent is given in figure 4. Large percentage errors are evident, the approxi- 
mate theory, in general, overestimating the lowest natural frequencies. The corre- 
sponding normal displacement modes given in figure 5 are found by using the numerical 
procedure and exhibit noticeable deviations from the assumed sinusoidal form of the 
approximate solution especially as n increases. Nevertheless, the approximate solution 
still exhibits the general shell behavior with respect to both the modes and frequencies. 
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Based on the comparisons made in this section, it is believed that the approximate 
theory gives reasonable estimates of frequencies and predicts trends adequately for 
values of r between -5 percent and 5 percent. All analyses in the remainder of this 
report are limited to shells contained within this range. 

EFFECTS OF MERIDIONAL CURVATURE 

By using the approximate analysis, the influence of meridional curvature is 
assessed by performing a parameter study for shells with differing meridional curva- 
tures as the length-radius and thickness-radius ratios are held constant. Results of addi- 
tional calculations based on membrane theory and pure bending theory are presented to 
clarify the relative roles of membrane action and bending action as they affect the vibra- 
tion of these shells. Although the results obtained are for a specific class of shell con- 
figurations, they may provide insight into the general behavior of more complex doubly 
curved shells of revolution. 

In calculating the effect of meridional curvature on shell vibration frequencies, the 
length- radius and thickness- radius ratios are held constant while the curvature ratio k£ 
is varied. Although the mass distribution varies slightly as k^ is changed, the effect 
of this variation on the natural frequencies is negligible if meridional shallowness is 
maintained. Thus, any change in natural frequencies can be attributed solely to the 
change in curvature of the shell and the resultant shell stiffness change. 

Effective stiffness of a shell is defined here to be that resistance per unit deforma- 
tion which the shell provides in a particular vibration mode. The natural frequencies of 
a shell structure are closely related to the effective stiffness of the structure in the sense 
that as the minimum effective stiffness increases, the lowest natural frequency increases. 
The positive curvature shells exhibit a strong stiffening character, the lowest natural 
frequencies increasing as the curvature of the meridian increases. For example, the 
lowest natural frequency of the positive Gaussian curvature shell in figure 2(a) is 
150 percent higher than that of the cylinder, whereas in figure 2(c), the lowest natural 
frequency is over 300 percent higher. On the other hand, negative curvature shells 
exhibit large decreases in frequencies as the circumferential mode number changes; 
these large decreases indicate an appreciable loss in effective stiffness. 

The lowest frequencies for each circumferential wave number occur at the simplest 
meridional mode m = 1 for positive Gaussian curvature shells and cylinders. However, 
this condition is not necessarily true for negative curvature shells. This trend can be 
seen in figure 4 where the lowest frequency occurs for m = 2 and higher meridional 
modes are associated with lowest frequencies for n > 3. This figure shows an inter- 
esting phenomenon in the spacing of the frequencies. As n increases, the spacing 
between frequencies associated with successive meridional wave numbers decreases. 
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This behavior suggests that experimental resolution of individual natural modes would 
be difficult to achieve at the higher n range. 


In order to interpret the results for negative Gaussian curvature shells shown in 
figure 2, a more extensive parameter study was conducted. The lowest natural fre- 
quencies and associated modes of a series of freely supported shells with a length- radius 
ratio of 3 (that is, s/R = 3) and with the rise-length ratio r varying from -5 percent to 
5 percent calculated with the approximate theory are presented in figure 6. The lowest 
natural frequencies for the shells of figures 2(a) and 2(b) are located in figure 6 by the 
vertical dashed lines. As the positive curvature increases, the lowest natural frequency 
and circumferential wave number increase monotonically. As the negative curvature 
increases from zero, the minimum frequency decreases and occurs in a progressively 
lower circumferential mode. For each mode, there is a distinct minimum in the varia- 
tion of frequency with curvature and hence, there are specific values of curvature for 
which large decreases in effective stiffness occur. The lowest natural frequency of 
cylinders and positive curvature shells occurs for m = 1. However, an m = 2 merid- 
ional mode is associated with the lowest natural frequency for shells within a small range 
of negative curvature. 

Figure 7 is a compilation of nine plots of the same type as that of figure 6. The 
length- radius ratio ranges from I to 10 and the radius -thickness ratio ranges from 100 
to 1000. The minimum frequencies, in general, decrease as the length increases and as 
the thickness decreases. The reductions in stiffness in the negative Gaussian curvature 
range are more prominent for the thinner shells, but the effect is noticeable to some 
degree for all the shells. The minimum frequency for each branch of the envelope in the 
negative curvature range occurs at the same curvature for a given length independent of 
the thickness of the shell. This result suggests that the curvatures at which the minimums 
occur are related to membrane action. On the other hand, the magnitude of the frequen- 
cies at the minimums is highly dependent on the thickness; thus, bending action must be a 
prominent factor in the effective stiffness of the negative curvature shell in the region of 
these minimum frequencies. 

MEMBRANE AND PURE BENDING ANALYSIS 

The results of the previous section imply that membrane behavior is closely related 
to the large reductions in effective stiffness observed for negative Gaussian curvature 
shells. The membrane equations associated with these shells have been solved in closed 
form for the membrane natural frequencies by assuming freely supported edge conditions. 
The membrane natural frequencies are given by equation (7) in terms of the nondimen- 
sional meridional curvature k^, modal wave length ratio j3, and the circumferential 
wave number n. The frequencies plotted in figure 8 for particular wavelength ratios for 
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unstressed shells are shown to be continuous linear functions of which decrease to 
zero as approaches the value given by equation (8). Since the total mass and mass 
distribution are essentially constant as k| is varied, this decrease in frequency must 
correspond to a decrease in effective membrane stiffness. Therefore, for a given modal 
wavelength ratio, there exists a negative Gaussian curvature membrane shell with a non- 
dimensional curvature k^ (given by eq. (8)) which vibrates without developing any effec- 
tive membrane resistance (that is, in an inextensional mode). As the meridional curva- 
ture increases negatively from this critical value, the membrane shell regains its 
stiffening characteristics. 

The lowest natural frequency envelope given in figure 6 is replotted in figure 9. 
Vertical dashed lines locate the values of r at which the membrane theory predicts 
inextensional behavior for the modes associated with the lowest natural frequencies. 

These values of r very nearly coincide with the rise ratios locating the minimum fre- 
quencies of the envelope. Hence, these reductions in lowest natural frequency occur at 
curvatures for which membrane theory predicts inextensional behavior. This behavior 
is discussed in reference 8 for static analyses where it is shown that pure bending states 
can occur in negative Gaussian curvature shells when in-plane displacement freedom (u 
or v) is allowed at the boundaries, a situation which exists for freely supported shells. 
This behavior can similarly be shown to be true for all the minimums in figure 7. The 
higher frequencies on either side of these minimums must be due to the recovery of 
membrane stiffness for that wavelength ratio as the curvature changes as is shown in 
figure 8. 

Since the calculations show that the effective membrane stiffness is negligible in 
negative Gaussian curvature shells for certain wavelength ratios, the frequency at these 
critical combinations of wavelength ratio and meridional curvature must be governed 
almost entirely by bending stiffness. Thus, the minimum natural frequencies found with 
a pure bending (inextensional) theory should yield values near the observed minimums. 

To verify this statement, the pure bending solution was obtained and is given in equa- 
tion (9). Horizontal dashed lines in figure 9 locate the pure bending frequencies for the 
modes occurring in the lowest natural frequency envelope. In each case, even though the 
pure bending results are independent of meridional curvature, the pure bending frequency 
closely approximates the minimum frequencies obtained by using the approximate theory 
for corresponding m,n modes. Hence, it may be concluded that the large reductions 
in frequency observed in certain negative curvature ranges are due to a loss in membrane 
action and that for these critical combinations of modal wavelength and meridional curva- 
ture, the shell is very nearly vibrating in a pure bending mode. Furthermore, these com- 
binations may be predicted from a simple membrane equation (eq. (8)). 
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It should be noted that every negative Gaussian curvature shell vibrates in a nearly 
pure bending mode for some particular wavelength ratio; however, if the m or n num- 
ber is large, there may be enough effective bending stiffness present to maintain a natural 
frequency well above the lowest natural frequency of the shell. 

CONCLUDING REMARKS 

An approximate set of equations of motion governing the vibration behavior of 
nearly cylindrical shells with shallow meridional curvature is developed and a closed- 
form solution is obtained for freely supported edge conditions. A comparison of the solu- 
tions from the approximate theory to those obtained from a more accurate theory shows 
that within a certain range of meridional curvature, the approximate theory yields 
acceptable solutions. 

The effect of the meridional curvature on the minimum natural frequencies of nearly 
cylindrical shells with either positive or negative Gaussian curvature is investigated. 

The positive Gaussian curvature shells have minimum natural frequencies well above 
those of corresponding cylindrical shells. The minimum natural frequencies of the nega- 
tive Gaussian curvature shells generally are below those of the corresponding cylinders 
and evidence wide variations in value, large reductions in magnitude occurring at certain 
critical values of curvature. The corresponding membrane and pure bending equations 
are also solved for the same edge conditions. Comparisons of the membrane, pure 
bending, and complete shell analyses show that the very low frequencies at the critical 
values of negative curvature result from a condition of essentially negligible effective 
membrane stiffness and, as a consequence, are dependent almost entirely on the shell 
bending stiffness. The membrane theory affords a simple method of determining the 
modal wavelength ratio and curvature combination at which this essentially inextensional 
behavior can exist for a negative Gaussian curvature shell, whereas the pure bending 
theory gives a good estimate of the magnitudes of the associated frequencies. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., January 21, 1969, 

124-08-06-11-23. 
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APPENDIX A 


DERIVATION OF APPROXIMATE SHELL EQUATIONS 

In this appendix, a set of linear homogeneous equations is developed which governs 
the dynamic behavior of doubly curved axisymmetric shells with constant shallow merid- 
ional curvature of either positive or negative Gaussian curvature. To achieve this, the 
nonlinear strain-displacement relations given by reference 6 are specialized for a 
cylinder and modified slightly by using the criterion of reference 7. The cylinder is then 
given a small axisymmetric deformation, and the strains caused by this initial geometric 
perturbation are removed to establish an unstrained state in the resultant doubly curved 
configuration. The nonlinear equilibrium equations corresponding to this unstrained state 
are developed and linearized. 

Nonlinear Strain Displacement Relations 

The first- approximation nonlinear strain- displacement relations developed in ref- 
erence 6 are based on the assumption of small strains and moderately small rotations. 
For a cylinder with radius R and coordinates and positive displacement directions as 
shown in figure 10, the resulting equations become 



where the symbol ~ is introduced to indicate quantities measured with respect to an 
unstressed cylindrical surface. 
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To permit the method of solution which follows, the term must be modified 

U g 

so that the term - — is eliminated. The quantity K ta can be rewritten as 

4R2 


where 


<>)i 


5^,9 ■’),{" 

contains only the linear terms of e 




2R 


(A2) 


£0* 


It is shown in reference 7 that 


linear terms of the type e^/R (where is any physical linear middle- surface 


strain) multiplied by a numerical factor of the order of unity may be added to the linear 
expressions for bending strains without introducing an error greater than that originally 
introduced through application of the Kir chhoff- Love hypothesis. Hence, lc ^ may be 
defined by 



(A3) 


In addition, only the Von Karman type of nonlinear terms are retained in the strain- 
displacement relations (that is, nonlinear terms composed of products of derivatives of 
the normal displacement w). The final form of the cylinder strain-displacement rela- 
tions reduces to 
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e e 
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(A4) 


A set of strain- displacement relations governing the behavior of a shell of revolu- 
tion with a shallow meridional curvature is developed in a fashion similar to that of ref- 
erences 4 and 5 by the introduction of an initial displacement w 0 (£) to the cylinder 
middle surface. Let 

w(£,0,t) = w(|,0,t) + w 0 U) (A 5) 
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where w is the displacement measured from the initially deformed surface along the 
normal to the cylindrical surface. Since the meridional curve is assumed to be shallow, 
it is sufficient as a first approximation to represent the initially deformed curve by 


w 0 (£) = - 



(A6) 


where s is the total meridional arc length and R| is the constant radius of curvature 
of the initially deformed meridional curve. (See, for example, ref. 8.) This representa- 
tion restricts the subsequent analysis to shells of the form shown in the sketches in fig- 
ure 1. As a consequence of equation (A6), if R| is large, the meridional displacement 
measured along a tangent to the initially deformed surface u and the circumferential 
displacement v can be written as 



v = v J 


With application of equations (A5), (A6), and (A7), the strains caused by the total deforma- 
tion can be written as 
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The strains caused by the initial deformation w Q are 




e 



(U)! 


€ ^0 = 0 
Kfc = - —• 

K s o ‘ 0 
“(So ' 0 


(A9) 


The initially deformed system is now taken to be in an unstrained equilibrium state. Any 
deformation u,v,w away from this new equilibrium state introduces strains given 
approximately by subtracting the initial deformation strain of equations (A9) from the 
total strains of equations (A8) ^for example, Hence, 


e 








6 ~ = I(v , u ,e, w ^ w ,e 

6 ?e 2V,« + R + R 


K 


r 







J 


(A10) 


where f is now measured along the meridian. 
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Equilibrium Equations 

Nonlin ear equations .- Equations (A10) are the strain-displacement relations for a 
shell of revolution with a shallow meridional curvature. To obtain a set of equilibrium 
equations consistent with equations (A10), the principle of minimum potential energy is 
applied. 

The potential energy n for an isotropic, linearly elastic, homogeneous, vibrating 
shell with a uniform lateral conservative surface load p is 


" ' C C C/2I?* + zj «) 2 + + z “ e ) 2 + 2> f « + Z ~ K $° + 


o] r%2 r ^ 

2(1 - V-)(?£q + R dz d0 d£ - p J J wR d0 d£ 


I ^ Ij l2 C Ca^- 1 ) 2 + + (*•*)! 


R dz d0 d£ 


(All) 


Performing the integration over z and requiring the variation of the energy to 
vanish yields 


611 = 



6e^ + N 0 6e e + 2N^ e 6e^ + M| 6 + Mg S/c^ + 2M £g 6 /c^ 


where 


- P 


6w - Kh^u t t + v t 


6v , + w t 

, L , L 



R d0 d£ = 0 


N^ = B(^ + Me~ e ) 

N 0 = B (e 0 + M€ ? ) 

■^£0 — 

) 

+ IXK^j 
M d = D^i< e+ iiK^ 
M^ e = D(1 - ii)k^ q 


(A 12) 


(A13) 
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with 


B = 


D = 


Eh 

1 - p2 

Eh 3 


~\ 


12(l - ul 2 ) 


(A14) 


Substitution of equations (A10) into equation (A12) and integration by parts in equa- 
tion (A12) leads to the following nonlinear equilibrium equations: 

Nt fl n ^ 
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£ £ £ “*■ 1 _ *• - — h ( N t + N £ n — *- 

R p2 R \ R 


w 




1 a w ,0 


+ rV N 0 -r" + N ^ w > V + p ■ ^ hW ,tt = 0 


(A15) 


,e 

with the following quantities prescribed at the boundaries where £ is constant: 

~\ 


N> 


2M 


ke 


"V'-R 

2Mfc fl a ~ ~ ~ w ft 

M> > + + NtW t + Nt fl ~>- 


M; 


or u 


or v 


or w 


> (A16) 


or w 




Linear equ at ions .- These equations are linearized for prestressed linear vibra- 
tions by the same procedure that was used in reference 3. Let 

N£ = n$($) + N 5 (£,e,t) 

Ne = Ng(4) + N 0 (£,e,t) 

^^0 = N |0^> 




(Equations continued on next page) 
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= M|(4) + M^,0,t) 
M 0 = M 0 (Q + 

Also let 

u = u(|) + u(|,0,t) 
v = v(|,0,t) S 

w = w(4) + w(i,e,t) 


(A17) 


(A18) 


where the barred terms are associated with the stresses and deformations due to an 
axisymmetric static loading and the unbarred terms are associated with stresses and 
deformations due to a subsequent infinitesimal linear vibration about the prestressed 
state. The equations governing the linear vibration about the prestressed state are found 
by introducing equations (A17) and (A 18) into equations (A15), applying equilibrium of the 
prestressed state, and neglecting nonlinear vibration terms. (See, for example, ref. 3.) 
The resulting equations are 

~\ 

Nt f + - vha ,. = 0 

R ,tt 


N 


6) 0 . XT , V , ~ 


R + _r2 

2M 


M 0,0 . 2M 1M _ 


R 


^ hv ,tt = 0 


M 




^^0 + ^ ^ + NtW 1 1 + —I w - ^hw .. = 0 

R R 2 R^ R ? r 2 ,00 ,tt 


(A 19) 


where the assumption = 0 has been made (that is, meridional static rotations are 

neglected) and where it is assumed that the static loading is such that = 0. These 
assumptions lead to constant membrane static stress resultants N| and Ng. 

The corresponding boundary conditions at £ = Constant are 


N fi = 0 


2M *0 

N ^ + -^=0 


R 


~\ 

or u = 0 
or v = 0 ( 


(Equations continued on next page) 
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2M 


M 


4,4 + 


£^ + n { w , 5 '° 


M£ = 0 


or w = Ol 
or w = 0 

J 


(A20) 


The equilibrium equations (eqs. (A19)) and boundary conditions (eqs. (A20)) apply 
to shells of revolution with shallow meridional curvature but have no such limitation in 
the circumferential direction. The distinction between these equations and a Donnell- 
type formulation (shallow in both directions) involves the underlined terms in equa- 
tions (A19) and (A20) which would be omitted for a Donnell-type approximation. The 
terms which would be omitted for the Donnell approximation are similarly underlined in 
all subsequent derivations. 


The vibratory stress_and moment resultants are found from equations (A10), (A13), 

(A17), and (A18). Since — = 0, the resultants are 

d4 


where 



+ / ie 0 


N e = B ( 6 e 
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6 R R 
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e 4e B 2VR 

> 
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e -) 
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(A21a) 


(A21b) 
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With the aid of equations (A21), equations (A19) are written in terms of displace- 
ments as 


i - m u ,ee 1 + m /J_ . iA 


v h 




M 2 r 2 2 R R/ " B ,tt 


r r 21 ja/ 2 r T rr,«« ir RifR 


* 2 (2 - m)a 2 £h v _ 0 

12 r 2 12 ,U0 B ,tt 


> (A22a) 


— + i±\u . + /i + _A 


* 2 (2 - M) 


X 2 4 

v + V w 


R^ + Ry^ + (r + rJr " !2r2 V ^00 ' 12 v ,^9 + 12R 2 


+ '^2 + ^ + 4 W -^ W ”-^ ! ^ + -»“ = 0 


lRt 2 T ^R T i2) w B B-^2-^ w ,tt 


J 


where 


X = 


R 


Similarly, the corresponding boundary conditions at | = Constant are 




u , + g- + ix ue + 2 . = o 

,1 Re R R J 


V + (‘♦Skr 


X_ 

12 


“(2 - + R 2 w + (2 - p)w 




N t 

•f,r 0 


or u = 0 


or v = 0 


or w = 0 


> (A22b) 


w ,«£ + ^( w ,o- i ), s = 0 


or w j. = 0 

1 S 


J 


The constant membrane prestress stress resultants are defined by 

= N 

N 


N 0 = R p - 


(A23) 


R 


D 
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where N is an applied meridional stress resultant and p is the applied constant 
lateral internal pressure. Equations (A22) are partial differential equations with constant 
coefficients. If prestress deformations had been retained, additional terms with variable 
coefficients would occur and numerical methods of solution similar to those used in ref- 
erence 3 would be required. 
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SOLUTIONS OF APPROXIMATE SHELL EQUATIONS 


General Solution 


A closed-form solution to equations (A22) is available for the freely supported 
boundary conditions (N^ = v = w = = 0). The solution can be expressed as 


u = U mn cos SHIJl sin nd e iu,t 


v = Vmn sin cos n0 e 1 


w = W mn sin 2^-=- sin n 6 e ia,t 


where m is the number of meridional half-waves, n is the number of circumferential 
waves, and u> is the natural circular frequency. This solution satisfies the freely sup- 
ported conditions at £ = 0 and £ = s. Equations (Bl) are substituted into equa- 
tions (A22) to yield a set of linear homogeneous equations. For a nontrivial solution to 
exist, the determinant of the coefficient matrix of the resultant set of equations must 
vanish. This procedure leads to the characteristic equation 

-ft 6 + A 2 ft 4 - A x ft 2 + A 0 = 0 (B2) 

where £2 is a nondimensional frequency parameter and where the coefficients Aq, Aj, 
and A 2 are as follows: 

A 0 = b llf 3 22 b 33 ■ ( b 2 3 ) " b 12( b 12 b 33 " b 13 b 23) + b 13( b 12 b 23 _ b 13 b 22) 

Ai = b 2 2 b 33 - (b 23 ) 2 + b 33 b 11 - (b 13 ) 2 + bi 1 b 22 - ( b 12 ) 2 ) ( B3 ) 


where 


A 2 - t>ii + b 22 + b 33 


v. /mirRr . 1 

b ll= l-i-J + “ 


b 12 = 


_ 1+jx rnwR 


bio = - 


'( k ^ + ^ 
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b 22 



+ 


g 


b 23 = -n (l + Mk|) - ^2 n 


n 2 + (2 - #i)( 


m7rR' 


h 

b 33 - 12 




+ k { 2 + 2Mk 5+ l + |(H|S) + ^„2 


R 

and where the ratio of curvatures has been introduced as kt = . The natural fre- 

quencies of a freely supported doubly curved shell of revolution with shallow meridional 
curvature can be found for any specific mode by solving equation (B2). This equation 
leads to three frequencies for each set of m,n considered, since in-plane inertia terms 
have been retained. 


Membrane Solution 


It is of value to inspect the extreme case of zero bending stiffness (D = 0). The 
vibratory behavior in this case is associated with only the extensional properties of the 
shell. The membrane equations are found directly from equations (A19) by deleting all 
moment terms to yield 


N 


M 


+ ^=o 

R 


n m 

R 




= 0 




N| Ne 

R| + R 



- N fl - ? ■- - + I'hw 
6 R2 



(B4) 


where, for convenience in this limiting case, the in-plane inertias have been neglected. 
The first two equilibrium equations are identically satisfied by the introduction of the 
stress function t// defined by 


r2 B 


B ~ R 



1-m 

R B 2 




J 


(B5) 
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The third equilibrium equation becomes 




, 96 


+ Ri/^ 


b 


N 


Raw ,S«-^‘ 


b w ,ee + ^ w ,tt = 0 


(B6) 


It can be seen from equation (B6) that for a negative Gaussian curvature negative) 
unstressed membrane shell, the governing equation has a hyperbolic character. The 
equation for the cylinder has a parabolic character and for the positive Gaussian curva- 
ture shell has an elliptic character as noted in reference 8 for static membrane theory. 

Equation (B6) is satisfied by 


4/ = 


membrane 


sin sin ne e iu>t 




w - ^membrane s * n s ^ s * n 


(B7) 


where it follows from equations (B5) that 


u “ ^membrane cos s ^ s * n e 


iwt 


v - V m embrane sin s ^ cos nd 


(B8) 


This solution satisfies the freely supported boundary conditions. 

With the application of equations (B7) and (B8) to equations (B5), equations (B5) 
may be written in matrix form as 


m7rR 
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+ /i 
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mTTR 
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2n( m7rR ) 
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membrane 


L 1 - M J 


where from Cramer's rule 


^membrane - ~ 
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{Zff + 

n 2 k| 


(B9) 


(BIO) 
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Substitution of equations (B7), (B8), and (BIO) into equation (B6) yields the following 
membrane frequency equation: 


^membrane _ 



(Bll) 


where /3 = ns is the ratio of the circumferential to the axial wavelength of the vibra- 
tion mode. The frequency determined from equation (Bll) vanishes for an initially 
unstressed membrane for meridional curvatures given by 

k| = -/3 2 (B12) 

Thus, it appears that for certain mode shapes the vibration of a negative curvature mem- 
brane shell with freely supported edges is not sustained by the membrane stiffness of 
the shell. This condition would suggest that the shell experiences inextensional behavior. 
The proof that this behavior occurs can be demonstrated by setting the middle surface 
strains of equations (A21b) equal to zero and substituting in equations (A21b) the solution 
given by equations (Bl). The resulting conditions for inextensibility are that 


= - 0 2 

^mn _ 1^ 
W m n' n > 

^mn _ _ j8 

Wmn " " n , 


(B13) 


It can be shown that these three conditions do occur when the membrane frequency is 
zero. 


Pure Bending Solution 

With the assumptions that B remains finite and that the middle surface strains 
equal zero, the pure bending (inextensional) equations are found from equations (A19) 
and (A21a). The resulting equations in terms of displacements can be written down 
directly from equations (A22a) by retaining all terms containing A 2 , and v and w 
inertia terms. Since in-plane displacements contribute mainly to in-plane stretching 
and, in general, have a negligible effect on inextensional deformation, a further simpli- 
fication can be made by deleting the v inertia and the non-Donnell terms (that is, by 
deleting the second equation of (A22a) and by deleting all derivatives of v in the third 
equation of (A22a)). The following equation results: 
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V 4 w + w , , = 0 (B14) 

The meridional curvature 1 /R| does not appear in equation (B14); thus in the approxi- 
mate formulation the pure bending state is independent of R^. 

The freely supported boundary conditions are satisfied by 

w = Wv, sin sin n 0 e^* 3 ^ (B15) 

o 


where the subscript b denotes bending. 

The natural frequencies of pure bending are given by 

V - til 1 + » 2 ) 2 " 4 

where 

R 2 w b 2 ^(l - M 2 ) 

V = — S 

and where again 


0 = 


ni77R 

ns 


(n * 0) (B16) 


These frequencies are equivalent to the natural frequencies of a simply supported rec- 
tangular plate with aspect ratio 7 tR/s. 
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(a) = ±0.05; s/R = 3; R/h = 1000. 

Figure 2.- Comparison of minimum frequencies from approximate analysis with those from numerical analysis for freely 

supported shells. 
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(b) = ±0.1; s/R = 3; R/h = 1000. 

Figure 2.- Continued. 
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± =1000 
h 

Numerical solution { ref. 3 ) 
° K *= .15 ( r = 5.6 2 % ) 

° K * = -/I5 ( r=~5.6 2 % ) 

— a — Approximate solution 



4 6 8 10 12 


Circumferential wave number , n 

(c) k£ = ±0.15; s/R ~ 3; R/h = 1000. 
Figure 2.- Concluded. 
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Figure 4.- Comparison of minimum frequencies from approximate analysis with those from numerical analysis for freely 

supported shells, = -0.05; s/R = 10; R/h = 1000. 
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Minimum nondimensional frequency , 
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Figure 8.- Effect of meridional curvature on the natural frequencies of freely supported unstressed membrane shells. 

k^ = -0.5 to 0.5; |i = 0.3. 






Minimum nondimensional frequency , 


CJ 



Percent meridional rise , t 


Figure 9.- Comparison of membrane and pure bending results with complete approximate theory, x = -5% to 5%; s/R = 3; 
R/h = 1000. The meridional wave number m is 1 unless otherwise noted; the circumferential wave number n is given 
on curves. All calculations based on approximate theory. 



National Aeronautics and Space Administration 
Washington, D. C. 20546 


POSTAGE AND FEES PAID 
NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION 


OFFICIAL BUSINESS 


FIRST CLASS MAIL 


O'! H 
i ' 

. i n 


ij 1 57 : >i 3 l)S 

. * ; C r *F A i-’Cls 3 
I ) r. I r H ^ e, ~ 


69086 00903 

L A'U it A r ( ■ R V / AFWL/ 
i6\S L ? * 1 1 W i* I X 1 - C 6 


7 11 f 


CHI! 


U i 


POSTMASTER: 


If Undeliverable ( Section 158 
Postal Manual) Do Not Return 


f( The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the tuidest practicable and appropriate dissemination 
of information concerning its activities and the results thereof” 

— National Aeronautics and Space Act of 1958 

NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and 
technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 
Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include conference proceedings, 
monographs, data compilations, handbooks, 
sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and Notes, 
and Technology Surveys. 


Details on the availability of these publications may be obtained from : 

SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 



